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INTRODUCTION AND SUMMARY 


This report discusses the two key items in the navigation signal 
subsystem of the NAVSTAR satellites, the L-band transmitter and ultra- 
stable oscillator or clock. The NAVSTAR interim report (Reference 1) 
proposed a 50-watt output solid-state design for the transmitter and a 
quartz crystal reference for the oscillator. This report discusses a 
refinement of that transmitter design along with two additional designs 
and a commercially available crystal oscillator which meets the 
NAVSTAR requirements and can be used for NAVSTAR satellite test 
programs . 

Of the two additional transmitter designs, one is a straightforward 
approach using an L-band TWT power amplifier and the other incorpo- 
rates a solid-state L-band power amplifier. This latter design promises 
to be the most attractive from all standpoints of design including relia- 
bility and DC-to-RF efficiency. Efficiency is expected to be greater than 
30 percent. However, the design also contains an element of uncertainty 
since development depends on the availability of high-power integrated 
modules or chips. Work is now in progress on high power S-band module 
at TRW Semiconductors. Present plans are that 20-watt modules will be 
developed in six to nine months; the design assumes the use of these 
modules scaled to operate at L-band. 

The first solid-state design (basically the same as in the interim 
report) using state-of-the-art components is shown to be superior to a 
TWT design in all respects except DC-to-RF efficiency. An efficiency 
of 22 percent minimum is expected as compared to 28 percent minimum 
for the TWT. 

A crystal oscillator is available from Frequency Electronics of 
New Hyde Park, N. Y. , which has shown the capability to meet the 
stability requirements needed for the NAVSTAR system. Similar units 
have been flown in space and extensive data is available which support 
the stabilities claimed for the unit. 



A test program is .proposed which will verify the performance of 
the Frequency Electronics unit. In particular, a set of phase stability 
measurements are proposed from which variate difference statistics can 
be derived. These statistics demonstrate the expected time error in 
extrapolating oscillator drift as proposed in the NAVSTAR system. 
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The design of the L-band transmitter was investigated with the 
goals of maximizing DC-to-RF efficiency and reliability. Two solid- 
state designs and a TWT design were studied. The first solid-state 
design is based on existing components and circuit techniques. The key 
parameters of this design are compared to a TWT design. An S-band' 
tube was assumed which can be easily scaled down to operate at L-band. 

The second solid-state design is based on a projected solid-state 
power device which should be available in six to nine months. This pro- 
jection is based on work being conducted at TRW Semiconductors. The 
basic difference between the two solid-state designs is the frequency at 
which the power amplification occurs. In the first approach, power 
amplification is at 500 MHz and a varactor diode tripler is used to in- 
crease the frequency to L-band. This approach is essentially the same 
as that proposed in Volume IV, Section 3.2 of Reference 1. In the second 
approach, power amplification is directly at L-band. This design has 
superior DC-to-RF efficiency because it avoids the loss in the varactor 
multiplier. 

Table 1 lists the design parameters for the transmitters. A nomi- 
nal carrier frequency of 1560 MHz is shown although the transmitter can 
operate anywhere in the 1540-1660 MHz band of interest. The 5 percent 
bandwidth specified is more than adequate to handle the signal modula- 
tion. The 50-watt output power is the minimum desired. Expansion to 
100 watts is considered as growth potential. 

The junction temperature of the transistors is restricted to 12 5°C 
maximum for reliability. This temperature is sufficiently degraded 
from maximum junction temperatures for silicon devices (200°C) to 
provide conservative design for reliability and is a generally accepted 
value for space application. 

Not shown in Table 1 but nevertheless a requirement on the trans- 
mitter is the time division transmissions from the satellites. That is, 
the transmitter must be gated "on" and "off" with control signals. During 
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Table i. Transmitter Design Parameters 


Carrier frequency 

1560 MHz 

Bandwidth 

5 percent or better 

Power output 

50 watts (growth to 100 watts) 

Maximum junction 
temperature 

125°C maximum 

Baseplate temperature 

50°C maximum 

Input power from driver 

1.6 watts at 520 MHz (design 1) 
63 milliwatts at 1560 MHz 
(design 2) 

50 milliwatts at 1560 MHz (TWT) 

Input signal modulation 

Phase modulation (PSK) 

Input and output impedance 

50 ohms 

Overall DC-to-RF 
efficiency 

Goal: 22 percent minimum 

(design 1) 

31 percent minimum 
(design 2) 

28 percent minimum 
(TWT) 

Power supply 

+28 V, +2 percent 

Input /output YSWR 

1. 5 maximum without breakup 

periods power must not be drawn from the power supply. This 


requirement is' easily met by interrupting the signal input to the solid- 
state transmitter and the beam voltage on the TWT transmitter. 

As a starting point m the design of the transmitter, an assessment 
must be made of the frequency at which power amplification should take 
place. The higher the frequency can be made the more efficient the 
transmitter should become because of the reduction in higher power level 
varactor multiplier stages. Amplification directly at L-band is preferred 
since high power varactor stages are eliminated. 

The best available solid-state power device at the present time is 
the 2N5178 of TRW Semiconductors rated at 70 watts at 500 MHz at 25°C. 
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Presently available transistors that are capable of power gain at L-band 
have poor thermal resistance (0.)* Consequently, the power output per 
transistor is limited to approximately 1 watt if the junction temperature 
is to be kept below 125°C, Therefore 50 transistors in parallel are 
required to produce 50 watts output. 

Under development at TRW Semiconductors are integrated circuit 
S-band power modules developing 10 watts. It is projected that in approxi 
mately six to nine months 20 watts per module can be achieved. 


2. 1 SOLID-STATE WITH VARACTOR MULTIPLICATION 

.The first solid state design utilizes the 2N5178 power transistor 
operating at 520 MHz and the Varian LSP 811A-3 varactor diode as a 
tripler stage for increasing the frequency output to 1560 MHz. The avail- 
able power output per transistor at 520 MHz is derived as follows. The 
thermal resistance 0. of the transistor derived from the transistor speci- 
fications is 


0 . 

J 


maximum junction temperature - ambient temperature specified 

power rating 


= 20 ° J 25 = 2. 5° C /watt 

The maximum baseplate temperature is 50 C. An additional 10 C tem- 
perature gradient between the transistor case and baseplate is allowed, 
so that 

Device temperature rise = (maximum allowable junction temperature) 

- (maximum baseplate temperature) 

- (allowable temperature gradient) 

= 125 - 50 - 10 = 65°C 


The total device dissipation is 


p _ total temperature rise 

d fT. 

3 


65 

2.5 


26 watts 
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The RF power output is given by 


where 



24. 3 watts 


G = stage gain = 2.5 (4 db) 

7 ] = collector efficiency = 60 percent 

Therefore, the total available output power per transistor at 520 MHz is 
24. 3 watts while maintaining the maximum junction temperature at 125°C 
with a baseplate temperature of 50°C and allowing an additional 10°C for 
the temperature gradient between the transistor case and the baseplate. 

The efficiency of the times three multiplier is given by: 


V 

and y 
where n 



harmonic number 



The best available varactor diode, the Varian LSP 811A-3, has 
f c = 90 GHz. Therefore, M = 0. 07 for 

f = 1560 MHz 
o 

f = 90 GHz 

c 

n = 3 

and the efficiency 7 ] is 86 percent for a loss of 0.7 db. The overall 
efficiency of the x3 multiplier, including the input, output and idler cir- 
cuits, is assumed to be 63 percent, equivalent to a 2 db loss in the 
multiplier. 
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The diode thermal resistivity is given as 4. 7°C per watt. Hence, 
maximum power dissipation per diode is 


maximum allowable temperature rise 
thermal resistance 


65 
4. 7 


- 13. 8 watts 


The minimum number of transistors and varactor diodes required 
for the transmitter are calculated below. 


Input power to 
x3 multiplier 

Number of diodes 


50 -watt Output 


100-watt Output 


50 

0763 


80 watts 


100 

0763 


159 watts 


{1 - 0.86) (80) 
13.8 


(1 - 0.86) (159) 
13.8 


1.6 


Number of 
transistors 


80 

24. 3 


3. 3 


159 
24. 3 


6 . 6 


With the aid of the preceding analysis, the following design approach 
was selected: 

1) 50 watts at 1560 MHz will be obtained by combining the 
output power from two identical paths, each delivering 
40 watts minimum at 520 MHz. The 40-watt outputs will 
be derived from three transistors in parallel. Two 
transistors would suffice; three are used for a conserva- 
tive design and to keep the junction temperature below 
125 C for increased reliability. 

2) A X3 multiplier containing two varactor diodes in each 
path delivers 25 watts minimum at 1560 MHz to the 
combiner. A balanced configuration is used in the multi- 
plier to reduce odd harmonics. Combining the two out- 
puts at 520 MHz and then multiplying was rejected 
because of increased design problems with the multiplier 
at the higher power levels. 

3) Growth to 100 watts at 1560 MHz is possible by using 
four or five transistors in parallel in each path delivering 
80 watts minimum at 520 MHz. This power level is not 
recommended because paralleling more transistors in- 
creases balancing problems in combining the outputs. 

Thermal design and arcing at the output become more 
severe problems especially under the constraints of 
packaging size and weight. 
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Figure i shows a simplified block diagram of the entire transmitter. 
The driver and phase modulator are straightforward. The driver delivers 
1.6 watts at 520 MHz to the solid-state power amplifier. The DC power 
input to the modulator and driver will be approximately five watts. 

Figure 2 is a detailed block diagram of the power amplifier showing the 
gain and power levels at each stage. This design uses a minimum number 
of active stages. 


1.6W 

I 

I 


MODULATION 
INPUT 



104 MHz 


Figure 1. Simplified Block Diagram of Transmitter 

The phase -modulated 520 MHz carrier from the driver is applied 
to the input stage Al, which is a single stage Class C amplifier exhibiting 
approximately 5 db of gain. The output feeds a 3 db power divider which, 
in turn, excites two identical transistor chains. A2, similar to Al, 
raises the power level to approximately 8 watts'. The power divider that ' 
follows excites a Class C doublet amplifier, the output of which is com- 
bined in the hybrid to produce over 18 watts. This, in turn, provides an 
input to a 1:3 power divider feeding the Class C triplet. The combined 
outputs in the hybrid that follows yields over 40 watts. The basic circuit 
for each transistor stage is shown in Figure 3. 

The X3 unit utilizes two varactor diodes in a balanced configura- 
tion. The two multiplier outputs at 27. 5 watts and 1560 MHz are then 
combined to produce 54 watts at L-band. The unwanted outputs are 20 db 
down at this point, but if greater suppression is required, a bandpass 
filter (three -pole will produce approximately 36 db of suppression with 
a bassband loss of 0. 3 db) can be added after the output combining 
network. 
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2.5 W 8.1 W 18.8 W 43.4 W 27.5 W 

520 MHz I I +4 db ' +4 db I I 



DC POWER INPUT = 243 W 
RF POWER OUTPUT = 54 W 
DC/RF EFF = 22% 


Figure 2. Block Diagram, Power Amplifier 












+28 V 



OUTPUT 
520 MHz 


The total number and types of transistors and varactor diodes re- 
quired are as follows: 


Stage 

Number 

Type 

A1 

i 

PT 6680 

A2 

1 X 2 

2N5177 

A3 

2X2 

2N5177 

A4 

3 X 2 

2N5178 


13 

Total transistors 

x3 

2 x 2 

LSP811A-3 


4 

Total varactor diodes 


The overall DC-to-RF efficiency of the transmitter is estimated 

to be not less than 22 percent; depending upon the availability of efficient 

high frequency varactor diodes, the design goal is 25 percent. The 

transmitter requires a +28 volt regulated (2 percent) power supply. 

Excluding output filter, the transmitter is estimated to weigh 4. 5 pounds 

3 

maximum and can be packaged in 12" X 5" x2" or 120 in . 

The following are the critical design areas for the transmitter. 
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1 ) 


A typical power combining network is shown in Figure 4. 
This technique allows a summation of "n" ports. The 
design is quite simple, yet care must be taken to insure 
identical path lengths for each input leg. To validate 
this design a power triplet (Figure i. Stage A4) has been 
designed, built (Figure 5) and tested to deliver 50 watts 
at 500 MHz. The transistors used were prototype units 
of the 2N5178. Test results are given in Table 2. In 
Figure 5 quarter wavelength coaxial transmission lines 
are coiled in order to save space. The main feature of 
this design is that transistors operating in parallel are 
isolated from each other by a hybrid. 

2) The x3. multiplier requires particular attention in order 
to find the optimum varactor diode to maximize the 
efficiency. The losses in the input, output, and idler 
circuits must also be minimized. The design is based 
on a 50-ohm input/output impedance so that the unit can 
be aligned and tested independently of the entire trans- 
mitter, diminishing the interface problems. 

3) Another area that will require particular attention is 
the final power combiner. The correct summation of 
the two multiplier outputs depends upon phase coherency 
of the signals. Final phase adjustment on one of the 
combiner inputs will be necessary. 


)/4 TRANSM LINE 
Zo = iJh 

N = NUMBER OF INPUTS 
2L 

*A\ — | 

OUTPUT 

INPUT 3 

3 INPUT COMBINER NETWORK 


INPUT I 
R = 

INPUT 2 

2 INPUT COMBINER 

Figure 4. Power Combining Networks 
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of Power Triplet 



Table 2. Test Results on Triplet Unit 


Transistor 

Number 

Power 

Input 

(watts) 

Power 
Output 
(watts ) 

I 28 V 
(amps ) 

Gain 

<db) 

Efficiency 

(percent) 

Individual Transistors 





i 

7. 8 

20 

1. o 

4. 1 

71.5 

2 

7. 7 

20 

0. 9 

4. 1 

79. 0 

3 

6.8 

20 

1. o 

4. 7 

71.5 

Triplet Unit 






i, 2, 3 

19.6 

50 

2.48 

4. 05 

71. 5 


2.2 SOLID - STATE WITH L-BAND POWER AMPLIFICATION 

The second solid-state design utilizes integrated circuit power 
modules operating directly at L-band and producing 20 watts maximum 
output power each. The modules do not exist at present but are pro- 
jected to become available by mid 1 9 6 9 • This projection is based upon 
development at TRW Semiconductors on 10 -watt integrated power chips 
operating at 2. 0 GHz. Three new devices which will deliver 1.0, 2.5 
and 5. 0 watts at 2. 0 GHz are now being introduced by TRW Semicon- 
ductors. Present plans are to emphasize development of increased 
power capability in the 1-2 GHz range rather than further extension of 
operating frequency at current power levels. 

Appendix A shows a projection made by the U.S. Army Electronics 
Command on development of solid-state devices. This projection was 
made in July 1968 and indicates the development of transistors with CW 
power levels exceeding 10 watts above 1 gHz by 19? 2* 

The concept for the transmitter is shown in Figure 6. The input 
from the phase modulator and driver is at 1560 MHz and requires only 
63 milliwatts. The power amplifier is composed of 10- watt and 20 -watt 
discrete integrated module stages and hybrid power dividing and com- 
bining units. The individual integrated amplifier modules utilize multi- 
transistor configurations on a single chip and micro striplines for internal 
matching and combining. Four 20-watt modules are shown for the output 
stage. The power per* module is derated to keep the junction temperature 
below 125°C. 
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6.3 MW 6.3 


1560 MHz, 


DRIVER 


MODULE 


•i-O uD 

rt>n 


20 W 
MODULE 


L-BAND 

OUTPUT 


(OPTIONAL) 


20 W MODULE 



TYPICAL 20W L-BAND INTEGRATED MODULE 


Figure 6. Design Concept for Direct L-Band 
Amplification Transmitter 


With microstripline within the module, considerable reduction in 
lead length is achieved, increasing the operating frequency and power. 
Up to 15 watts at 2 GHz has been demonstrated using these techniques. 
The individual modules are designed with a 50-ohm input and output 
impedance so that their performance can be checked prior to integration 
into the system. The hybrid combiner and divider networks can also be 
constructed utilizing the stripline approach. Use of the stripline has an 
important advantage in that later mistuning of the networks is unlikely 
since distances once established will remain fixed. 








This design seems more adaptable to 100 watts. The power can 
be extended by adding more modules in parallel. Matching and combining 
added units is not as serious a problem as for the first design because of 
the 50-ohm resistive input-output of each module and the microstripline. 

The DC-to-RF efficiency for the 50-watt unit is estimated as 
follows. The DC power input to the power amplifier is 


dc 


i + —L. + i 

G 1 G 1 G 2 


where 

= E .F power output 

y = collector efficiency of each stage 
Gj = gain of final stage (4 db) 

= gain of next stage (5 db) 

.The collector efficiency is conservatively estimated to be 46 percent. 
Overall efficiency of the power amplifier 


V 



i + 



1 


G 1 G 2 


X 100 percent 


=31 percent 


This efficiency depends primarily upon the power gain of the last stage. 
Therefore, particular attention will have to be paid to this area. 

It is difficult to assess the size and weight of the power amplifier 
at this point, although it is anticipated that both will be significantly 
lower than possible with a discrete component design. As an example, 
a 5-watt module at 2. 0 GHz measures 1. 5" x 1. 0" x 0. 5". 

Thermal resistance of these units is not known at this time, but 
good heat transfer from the active devices to the baseplate will be 
essential; hence a solid conducting mass for the housing will be required. 
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The following are the critical design areas for the power amplifier: 


1) Since the present development on the high power 
modules is at S-band, scaling down will be required 

2) Good familiarity with micro stripline techniques will 
be necessary so that size and weight reduction can 
be achieved 

3) Since maximum overall efficiency is a prime goal, 
considerable attention will be paid to the last power 
stage to optimize its performance 

4) It is anticipated that a substantial effort will be spent 
in the area of thermal conductivity and interfacing to 
keep the junction temperature below 125 C. 

2. 3 TWT DESIGN 

The third design of the transmitter uses a traveling wave tube 
operating at E-band to obtain 50 or 100 watts output. A phase modulator 
and driver nearly identical to that used with the second solid-state design 
is required. The driver output power is approximately 50 milliwatts for 
the 50-watt output, since a TWT typically has about 30 db gain. 

The significant problem with this design is the availability of a 
space -qualified TWT for E-band application. Also the TWT requires an 
associated power supply or converter which generates all of the voltages 
to operate the tube from the primary +28V source. 

Although considerable success has been obtained in the use of TWT 
transmitters for space applications, no TWT's at 50 or 100 watts at E- 
band have been flown. 

Approximately two years ago, both Watkins -Johnson and Eimac 
(a division of Varian) embarked on a JPE-funded program to develop a 
50 -watt S-band TWT. The program is completed, and Eimac now claims 
to have a 50-watt S-band unit, and Watkins- Johnson a unit capable of both 
50- and 100-watts at S-band. Although the units have wide bandwidths 
(nearly one octave) a “scaling down" redesign is necessary for E-band 
operation. 

Both of the companies estimate that the scaling down would cost 
approximately #100-150K (6-9 months redesign time) which would be a 
nonrecurring cost. In addition, a space qualification cost would have to 
be added, approximately $100K, depending upon the requirements. 
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The units are designed specifically for space application with an 
intended lifetime of 50, 000 continuous hours. Life expectancy at 90 per 
cent confidence level is estimated to be 20, 000 hours. The size and 
weight of the redesigned S-band units would increase approximately by 
an amount proportional to the frequency difference. 

2. 4 COMPARISON OF SOLID-STATE AND TWT DESIGNS 

The key parameters of the TWT and the first solid-state design 
are compared in Table 3. Realistic comparison with the second solid- 
state design are not possible because of the unknowns in the proposed 
devices. Preliminary indications of size, weight, and volume compare 
favorably with those shown in the table, and an expected DC-to-RF 
efficiency of at least 30 percent makes the second solid-state design 
attractive. Moreover, the use of integrated circuits and microstripline 
construction should result in its being the most reliable of the three. 

The MTBF in the table are calculated in Appendix B. The calcu- 
lations indicate that the solid-state transmitter using normal failure 
rates for the components has a higher MTBF than the TWT transmitter. 
It is perhaps worthwhile pointing out that the capacitor failure rate used 
in the solid state calculations occurred during alignment and test but not 
during orbital operations. 

Table 3. Comparison of Transmitter Designs 


Parameter 

Solid-State 
(Design One) 

TWT and Power Suppl - 
(Converter) 

Weight (pounds maximum) 

6. 0 

8. 5 

Size (inches maximum) 

12x5x2 

15x6x4 

Volume (cubic inches) 

120 

360 

DC to RF efficiency (% minimum) 

22 

28 

Bandwidth (%) 

£10 

M0 

Gain in db (maximum) 

As required 

30 nominal 

Power to weight ratio 

8. 3 

5.9 

Temperature range (°C) 

-20 to +65 

-30 to +100 

Vibration (operating) 


20 g's, 20-2000 Hz 

Shock (nonoperating) 


200 g’s, 1 msec 

MTBF (hours) 

39, 200 

32, 400 
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A figure of merit based on the following formula: 


Figure of Merit 


DC/RF eff. x MTBF 
weight x volume 


shows an advantage of 4:1 for the solid-state design. This comparison is 
not totally fair since the parameters in the figure of merit are not of 
equal importance. For example, the overriding consideration in a choice 
may be the DC-to-RF efficiency since the satellite is power limited. 
However, the figure of merit does suggest the attractiveness of the solid- 
state design as compared to the TWT. 
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3. ULTRA- STABLE SATELLITE OSCILLATOR 


It lias been determined that a crystal reference oscillator will be 
adequate for the clock reference on each of the NAVSTAR satellites (see 
Appendix C). A crystal oscillator has physical characteristics which 
make it the most attractive unit for satellite application. In addition, the 
high initial cost and questionable reliability of atomic standards make the 
crystal oscillator even more attractive. Table 4 compares some of the 
characteristics of different oscillators. Although the figure of merit in 
the last column is somewhat arbitrary, nevertheless it suggests the 
attractiveness of the crystal oscillator. 

Table 4. Comparison of Oscillator Characteristics 



Volume 

(ft 3 ) 

Weight 

(lbs) 

Power 
Cons umpti on 
(watts ) 

Stability 

per 

12 Hours 

Figure of Merit 
(product of the 
four) 

Atomic 

Hydrogen 

16.4 

800 

200 

2 x 10" 14 

2. 56 x 10 -8 

Rubidium 
Gas Cell 

0. 5 

38 

42 

5 x 10“ 12 

-9 

4. 0 x 10 

Cesium 

Atomic 

Beam 

1. 5 

64 

60 

2 x 10" 13 

1.1x10 ^ 

Crystal 

Oscillator 

0. 082 

6 

2. 5 

1 x 10 -11 

1.2 x 10“ U 


A survey of companies that can produce the high precision crystal 
oscillators required for NAVSTAR indicates that few companies have built 
such units. Fortunately, one company has built an oscillator which 
appears to fit the NAVSTAR requirements and which has flown in space. 

A specification largely based on this oscillator has been written; a test 
program is recommended to confirm the performance of the oscillator. 


3.1 CRYSTAL AGING 


One of the main problems with crystal oscillators is the settling 
time or frequency stabilization when the oscillator is turned on. In 
general, frequency stabilization of a high quality oscillator follows a 
predictable curve.- This stabilization is referred to as aging. However, 
it usually takes several months for a conventional oscillator to reach a 
stable frequency (Reference 2). Frequency stabilization with time follows 
a logarithmic law 

X (t) = a + blog e |\ + 

where a, b, and t are independent constants of the crystal. The aging 
rate is obtained by differentiating this expression 

b 

d /Af 
dt \^f 

o 

Therefore as t >>t Q the crystal frequency tends to stabilize. To mini- 
mize the settling time of the crystal t should be as short as possible. 

Aging is considered to result primarily from mass loading changes 
on the quartz plate by the absorption of the contaminants that are trapped 
during the crystal mounting. Reduction in contaminants will decrease 
both the settling time and the magnitude of the aging. 

Until recently, crystals have been encapsulated at normal tempera- 
tures, resulting in long stabilization periods. In addition, after equili- 
brium is reached and operation of the oscillator is interrupted, the aging 
time again follows the original characteristic (see Figure 7). 

A new method utilizing a high temperature (400°C) hydrogen envi- 
ronment encapsulation has been developed to bake out the impurities. 

This reduces the contamination considerably, and reduces the aging from 
several months to a few weeks (see Figure 7). In addition, the re- 
establishment time of stable frequency after disruption is considerably 
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Figure 7. Aging Characteristics of Quartz Crystals 




improved. The high temperature process is possible due to introduction 
of thermocompression bonding of the ribbons (the electrical contacts) to 
the quartz plate. Bake-out of solder-bonded units is limited to approxi- 
mately 125°C. 

With the present day crystal technology, a 5 MHz perpendicular 
field, fifth overtone crystal employing ''Hold-seal" and thermocompression 
bonding is the best unit in terms of frequency and temperature stability. 
Also it is the least susceptible to mechanical shock and vibration. 

3. 2 CRYSTAL OSCILLATORS FOR NAVSTAR 

Out of the companies surveyed involved in the design of high quality 
crystal oscillators, two were singled out on the basis of having demon- 
strated actual units flown in satellites: Applied Physics Laboratory, 

John Hopkins University, Maryland; and Frequency Electronics, Inc. , 

New York. 

The unit developed by the Applied Physics Laboratory was for the 
Dodge satellite. Its characteristics are as follows: 


Operating frequency 5 MHz, fifth overtone horizontal 

field crystal supplied by Bliley 
Electric Company, utilizing 
conventional glass encapsulation 


Temperature control 
Aging rate 


Drift rate (rms) 


Buss voltage sensitivity 


Double proportional oven 

2x 10“ 10 per day after 6 months 
of operation; 5 x 10“ after 
3 weeks of operation 

2 x 10” per day 
2 x 10~: per 2000 second 

1 x 10' 10 per 200 second 

2x 1 0 ” 1 1 / volt 


Power consumption 

Ambient temperature 
sensitivity 

Load resistance 
sensitivity 


1 . 5 watts 
2 x 10" U /°C 

1 x 10 ’^/ohm at 50 ohms 
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Load capacity 2x10 **/ pf at 60 pf 

sensitivity 

Temperature range -23 to +47 °C 

Weight 2. 8 pounds 

Dimensions 3. 6 inches diameter x 5* 0 inches 

high. 

The above unit is not commercially available since it was designed 
ander a Navy contract for a specific application. Actual in-flight data is 
rot available. However, it is known that one of the units had a reported 
thermostat failure after four months of orbital operation. No remote 
frequency adjustment is provided. The price would be approximately 
S7-9K plus space qualification. The tests and delivery time would be 
approximately 6-8 months. 

The unit designed by Frequency Electronics (Figures 8 and 9) was 
also developed for the Navy under USNRL Contract N00173-66-C-0162 
md ONR Contract N0014-66-C-5220. The following are the unit 1 s 
characteristics: 

Operating frequency 5 MHz, fifth overtone horizontal 

field crystal, utilizing high 
temperature thermocompression 
bonding and M Kold-seal" encapsu- 
lation (the crystal is manufactured 
by the company) 

Temperature control Triple proportional oven 

Aging rate 1-2x10 ^per day after 7 days 

of operation 

i-2 x p er 

1-2x10 ^ P er hour 
4-8 x 10 per second 

1 x 10 ^ for 5 percent voltage 
change 

Power consumption 2.5 watts, including remote 

frequency control 

Ambient temperature 1x10 ^/°C 


Drift rate 


Buss voltage sensitivity 
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NOT REPRODUCIBL 






NOT REPRODUCIBL 



Load resistance 
sensitivity 


i x 10 ^ for 20 percent load 
change 


Remote frequency 
sensitivity 


Temperature range 


Utilizes mechanically variable 
capacitors, fine and coarse, 
adjustment by gear control unit 


Fine control resolution +1x1 
Coarse control range + 2x 10“ 
Control signal “digital 


“11 


Present design of -20 to +45°C 
can be extended 


W eight 


4 pounds, 6. 5 pounds with remote 
tuning 


Dimensions 


3* 5 x 4. 5 x 9. 25 inches (with 
remote frequency control unit). 


The design represents the state of the art for space -qualified high 
stability oscillators. The unit is commercially available, and the vendor 
is willing to cooperate in conducting additional tests. The price for an 
actual deliverable unit is shown below. 


In quantities 
of two 


[ Without remote frequency adjustment 
With remote frequency adjustment 
Space Qualification tests 


$ 8. OK each 
$511. 5K each 
$12. OK 


An oscillator built by Frequency Electronics for the Naval Air 
Systems Command was flown on an experimental satellite early in 1967. 
The 5 MHz oscillator frequency was counted down to 100 kHz and then 

used to modulate the 400 MHz carrier. The demodulated frequency was 
then compared to a 100 kHz hydrogen reference at the ground station. 

The oscillator which used a double proportional oven and a Bliley crystal, 
exhibited a temperature coefficient of approximately 1. 7 x 10 ^/°C, 
which since then has been improved to a 1 x 10 ^/°C by using a triple 
proportional oven. Figure 10 shows data obtained on the oscillator drift 
rate with time. The discontinuity at approximately 350 days is un- 
explained. Since this unit was flown considerable improvement in initial 
drift rate (aging) and vibration sensitivity has been made by improved 
crystal preparation methods. 
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Figure 10, Oscillator Drift Rate Versus Time 
3. 3 RECOMMENDATIONS 

The Frequency Electronics unit is recommended for use in a test 
program to evaluate the NAVSTAR application. It appears that the 
NAVSTAR requirements (see Table 5) can be met based upon data obtained 
from orbital flights and ground tests. The vendor claims that the final 
deliverable units can be reduced somewhat in size, and the temperature 
range can be extended. Presently the crystal is operating at an inflection 
point of approximately +60°C, where the aging rate is minimum. Operating 
the crystal at a higher temperature can be done at the expense of increased 
long-term drift. 

The following tests should be performed on the oscillator: 
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1) General performance tests 

a) Output power 

b) DC power consumption 

c) Heater stability and thermal time constant 

2) Output waveform evaluation using harmonic and spectrum 
analyses 

a) Harmonic content . . 

b) FM index using a multiplier in order to increase resolut 


3) Stability measurements 

a) Phase stability - with an atomic standard reference using 
a high accuracy phase meter (measure time e jror a g ains 
accumulated time with a resolution of at least 0. 02pseconds) 
Frequency stability - by the use of multiplying chains an 


b) 


^ J J ■> 

7* atomic sta.nda.xd reference 


4) Remote frequency correction evaluation including accuracy 
response and settling time 

5) Environmental tests using the methods in 3 above 

a) Temperature 

b) Shock and vibration 

c) Sensitivity to power supplies and load impedance 

d) Humidity and pressure 

e) Radiation. 


The phase stability measurements are of particular importance to 
obtain the proper oscillator error statistics for the NAVSTAR system. 
From these measurements "variate differences" can be derived, which 
are used to determine prediction error in the oscillator drift. See 
Appendix C for details on the variate difference measurements and their 

significance. 
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Table 5. NAVSTAR Specification for the Crystal Oscillator 


Operating frequency 
Frequency stability 
Aging rate 

Temperature sensitivity 

Buss voltage sensitivity 

Power supply 

Power consumption 

Ambient temperature 

Manual adjustment 

Remote tuning (control signal to be 
specified later) 

Spurious outputs 

Harmonics output 

Vibration (exact specification to be 
supplied later) 

Dimensions 

Weight 

EMI 

Humidity and pressure 


5 MHz 

-12 

1 x 10 per 3 hours 

1 x 10 ^ per day after 10 days 
of operation 

i x 10“ i2 / P C or better 

1 x 10 / 1 percent voltage change 


+12V using 1 percent regulation 
2, 5 watts maximum 


-10 to +55 C 


,-12 


Stability of 1 x 10 

Fine tuning setability x 10 ^ 
Coarse range +3 x 10~^ 


Down not less than 80 db 


Down not less than 40 db 


The unit is to exhibit a stability of 
1 x 10~® during the launch period 

4. 0x5. Ox 10 inches maximum 


6. 5 pounds maximum, including 
remote tuning unit 

To be specified 

To be specified 
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APPENDIX A 

PROJECTION OF SOLID-STATE DEVICE 
DEVELOPMENT 


by 

U.S. Army Electronics Command 
July 1968 
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U. S. ARMY ELECTRONICS COMMAND 

Fori Monmouth, New Jersey 



ELECTRONIC COMPONENTS LABORATORY 


TECHNICAL GOALS AND PROJECTIONS 


I JULY 1968 



TECHNICAL COALS AMD PROJECTIONS HAS BEEN PREPARED BY THE ELECTRONIC COMPONENTS LABORATORY 
AND IS PUBLISHED FOR THE INFORMATION AND GUIDANCE OF ALL CONCERNED. THE PROJECTIONS AND COALS 
SHOW! ARE WHAT WE FEEL CAN PE ACCOMPLISHED IN VARIOUS MISSION RELATED SPECIALTY FIELDS IN THE 
> TIMS FRAMES INDICATED CONTINGENT UPON THE AVAILABILITY OF RESOURCES AND/OR PRIORITIES, 

w 



SOLID STATE DEVICES • 


This effort provides for the continuous exploratory development of solid state devices, including 
investigations of those materials and processing tecnniques which will lead to the development of 
advanced semiconductors or semiconductor-like devices . 

Transistors & Diodes 


> 
» i 


The goals to be achieved in the discrete semiconductor device area include higher output powers 
from high frequency transistors, i.c., 2CD to 500 watts PEP at 30 MKz, 50 to 100 watts at 500 MHz, 
and 25 watts at 1000 KHz. Three types of transistors are required; CW for Ft-I operation, linear for 
single sideband use, and pulse for radar and pulse type communications . In the' lower frequency 
semiconductor area audio units able to dissipate up to several kilowatts are required. In addition 
switching devices, which can switch several thousand amperes, will be developed for low voltage dc 
to dc convertors and pulse modulators. High input impedance, linear, small-signal field effect 
transistors for front end receiver amplifier use ranging to microwave frequencies; and large signal, 
power (tens of watts) lower frequency field effect transistors will be developed. 

Studies and ivestigations of improved semiconductor techniques, device improvement, new packaging 
concepts, and high power combining transistor circuit techniques will be advanced. They will include 
devices having higher breakdown voltages for AM applications in the UH? range: development of linear 
and varactor type transistors for SS3 and frequency multiplication applications; development of long 
thermal time constant power transistor for pulse power applications, thereby hardening the device 
thermally to applied transients; investigations of laminated silicon transistors utilizing their 
smaller devi.ee geometries, and advantages of higher power capabilities, higher frequency responses, 
and higher breakdown voltages; the investigation' of rate effects of four (4) layer structures; and 
low cost, low reactance packaging techniques, and evaluation. 


The investigation of noise, its cause and theoretical behavior in now semiconductor devices 
applicable to military equipments, will be continued. In addition the internal work will continue 
with the goal of characterizing, evaluating and providing assistance in the application of new non- 
microwave, non-mi croelectronic semiconductor devices to military equipment. On the second graphic 
"TRANSISTORS - Power Output vs frequency" the current curve indicates performance wherein power 
and frequency are limited, to a degree, by the material and processing methods. Future progress, 
i.c., through 1972 , is based on the premise of an advanced state-of-the-art vhere optimum performance 
is obtained for silicon and new and more sophisticated processing techniques are utilized. 



Unforeseen technological "breakthroughs and nev materials would make the prediction more optimistic. 

The third chart also covers transistors. The power gains predicted for CW and SSB transistors 
arc based on realistic "trade-offs'' involving both frequency and efficiency; that is, increased 
power gain would result if lower efficiency were acceptable. 

The field effect transistors, being essentially linear devices, will replace power vacuum 
tubes in 33B equipments at frequencies up to ?6 MHz. Without unforeseen breakthroughs, higher 
frequency operation or powers in excess of 60 watts are not likely before FT-72. 


> 

I 

LTt 



SEMICONDUCTOR DEVICES 


c 


ADVANCE OF TECHNOLOGY VIA APPLIED RESEARCH 


] 


RE LIAB I LITY 

IMPROVED METAL1ZATION TECHNIQUES 
SURFACE STABILITY INVESTIGATIONS 
ADVANCED DIFFUSION EPITAXY TECH. 
NEW MATERIALS 
NEW FABRICATING TECHNOLOGY 
IMPROVED THERMAL CAPABILITIES 
INFRARED EVALUATIONS 


I 


TRANSIENT 

INVESTIGATIONS 

HIGHER BREAKDOWN VOLTAGES 
STUDY OF MICRO 6e MESO PLASMAS 
HIGH TEMP. CONTACT MATERIALS 
PULSE THERMAL RESISTANCE 
TECHNIQUES 


I 


HIGH POWER - 
HIGH FREQUENCY 

NEW LEAD TECHNIQUES 
PROCESS IMPROVEMENTS 
IMPROVED PACKAGING TECHNIQUES 
ADVANCED SEMICONDUCTOR 
TECHNOLOGY 
ULTRA-RELIABILITY 


I 


l 1 — 

DESIRED SOLID STATE DESIGN CAPABILITIES 

p 


HIGHER EFFICIENCIES HIGHER FREQUENCY 

GREATER LINEARITY HIGHER MTSF 

INCREASED BANDWIDTH LOWER COST 

HIGHER FOWCR LOWER NOISE 



A 


cn 


IB 


HIGH POWER VHP* UHF SOLID STATE 
PROTECTIVE TECHNIQUES 
MATCHING TECHNIQUES 


USER 


VEHICULAR. FIXED tk ROTARY WING AIRCRAFT COMM. 
MANPACK. COUNTERMEASURES. SURVEILLANCE. RADAR 






ENVELOPE 'POWER ' POWER GAIN 

WATTS 


TRANSISTORS 


DEVELOPED OBJECTIVE 

... -J a a =3 t=3 a t=a 





OUTPUT POWER (WATTS) 



TRANSISTOR FAILURE RATE VS OUTPUT POWER 



FAILURE RATE - PERCENT / IOOO HRS. 




NOISE FIGURE (DB) 


TRANSISTOR NOISE FIGURE VS FREQUENCY 



FREQUENCY (GHz) 


APPENDIX B 


RELIABILITY CALCULATION FOR 50 -WATT TWT 
AND SOLID STATE TRANSMITTER 

Reliability predictions have been prepared for the 50-watt solid 
state transmitter (first design) to a 50-watt TWT and its converter. The 
results are presented in Table B-l. Table B-2 lists the failure rates for 
each part used in the predictions. Most of these failure rates are based 
on orbital operating hours accumulated by the OGO, Vela, and Pioneer 
Spacecraft. The low failure rate for power transistors is the maximum 
likelihood estimator based on the number of failures and part hours. The 
high failure rate for power transistors is the upper one-sided confidence 
limit at 60 percent confidence level. The low failure rate for variable 
capacitors is based on engineering judgement since it is a low usage part 
and has not accumulated many orbital hours. The high failure rate for 
variable capacitors is the upper one-sided confidence limit at 60 percent 
confidence level. The failure rate for the 50-watt TWT was supplied by 
the Eimac Division of Varian Associates. TRW has no orbital data on 
TWT’s of this power rating. 

All failure rates were adjusted by a computer program (TESRX) for 
the principal electrical stress and temperature on each particular part. 
The program adjusts the failure rates by means of a four-parameter 
Eyring model which has been fitted to the temperature-stress curves of 
MIL-HDBK-217A for each generic part type. 

To predict the reliability of a 50-watt TWT and converter, it was 
necessary to obtain a parts count for a typical TWT converter. The TWT 
converter in the Pioneer spacecraft was selected. 

To provide a range of values, three predictions were prepared using 
three different part case temperatures, 50°C, 60°C, and 70°C. For each 
prediction, the principal electrical stress on the parts was assumed to be 
30 percent on semiconductors, 50 percent on solid tantalum capacitors, 

25 percent on the other types of capacitors, and 25 percent on resistors. 
These values were selected as typical stress levels based on other TRW 
designed converters. The parts count and failure rates are listed in 
Table B-3. The results indicate a range from 28, 600 to 32, 400 hours 
MTBF for the TWT and converter. 
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Table B-i. Predicted MTBF's 


Failures MTBF = Mean Time 

Per 10^ -hours Between Failures {hours 


TWT and Converter 


At 50°C Case Temperature 

30, 887 

32, 400 

At 60°C Case Temperature 

32, 793 

30, 500 

At 70°C Case Temperature 

35, 013 

28, 600 

Solid State Transmitter 



With low failure rates for variable 
capacitors and rf power 
transistors. 

20, 129 

49, 700 

With low failure rate for variable 
capacitors and high failure rate 
for rf power transistors. 

25, 524 

39, 200 

With high failure rate for variable 
capacitors and low failure rate for 
rf power transistors. 

46, 979 

21, 300 

With high failure rates for variable 

5-2, 374 

19, 100 


capicators and rf power 
transistors. 
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Table B-2. Part Failure Rates 


Part Type 

* 

Failures/! 0^ Hours at 25% 
Stress and 30°C 

Diodes 

Varactor 

276 

Silicon Power 

14 

Zener 

60 

Transistor 

Silicon Power 

560 863 

Resistor 

Carbon Composition 

7 

Metal Film 

i . 4 

Wi rewound 

65 

Capacitor 

Ceramic 

15 

Mica, dipped 

3 

Glass 

23 

Tantalum Solid 

14 

Variable 

40 219 

Transformer 

Bow Voltage 

15 

High Voltage 

30 

Inductor 

15 

Traveling Wave Tube (TWT), 

20, 000 

50 watt 
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Table B-3. TWT and Converter 


W 

I 


Part Type 


Diodes, silicon power 
Diodes, zener 

Capacitors, solid tantalum 
Capacitors, ceramic 
Capacitors, mica 
Capacitors, glass 

Resistors, carbon comp. 
Resistors, metal film 
Resistors, wirewound 

Transistors, silicon power 

Inductors 

Transformers, low voltage 
Transformers, high voltage 

Traveling Wave Tube (TWT) 

TOTAL 


Quantity 


Total Stressed Failure Rate 
1 01 


Fai 


At 50 C Case 
Temperature 


a 


Hours 


At 60 C Case 
Temperature 


At 70 C Case 
Temperature 


48 

3 


1008 

270 


1152 

312 


1344 

360 


13 

7 

1 

4 


667 

168 

2 

192 


867 

238 

3 

296 


1022 

336 

3 

448 


15 

11 

1 


300 

22 

69 


540 

22 

74 


930 

22 

79 


10 

5 

3 

3 

1 


7970 

75 

45 

90 

20 , 000 

30, 878 


9070 

75 

45 

90 

20, 000 

32, 784 


10, 250 

75 

45 

90 

20 , 000 

35, 004 


Note: The following principal electrical stress levels were assumed for this prediction: 

Semiconductors — 30% 

Capacitors, solid tantalum — 50% 

Capacitors, all other types — 25% 

Resistors — 25% 

The case temperatures of the parts will depend on the baseplate temperature and the packaging 
technique. A tightly packed unit will cause higher case temperatures since less space is available 
for heat dissipation and vice versa. 



The block diagram for the solid-state amplifier is shown in Fig- 
ure 2 and the schematics of the amplifier stages and hybrids are shown 
in Figures 3 and 4. 

For this prediction, a baseplate temperature of 50°C was specified. 
The electrical stresses on the semiconductors were calculated using the 
power levels given in Figure 2. An electrical stress of 25 percent was 
assumed for all other parts. Four predictions were prepared by using 
the two different failure rates for power transistors and variable capa- 
citors discussed above. 

The power transistors proposed for the transmitter are new high 
frequency, high power devices. There is no data available to indicate 
the reliability of these devices. Therefore the low and high failure rates 
were selected to provide a range of predictions. The low failure rate 
would apply to transistors procured with a complete parts control program 
and the high rate would apply if there was a limited program. In either 
case the transistors would require PT specifications and qualification 
testing. 

The variable capacitors have not accumulated many orbital hours 
and therefore have a high failure rate at the 60 percent confidence level. 
The low failure rate is typical of what has been used in reliability pre- 
dictions for other programs at TRW. 

The parts count and stressed failure rates are listed in Table B-4. 
The results indicate a range from 19, 000 to 49, 700 hours MTBF for the 
solid-state transmitter. 

The solid state transmitter can be made more reliable by reducing 
the number of variable capacitors, replacing some of the variable capa- 
citors with fixed capacitors, or controlling the reliability level of the 
variable capacitors by instituting a complete parts control program. 

The tuning procedure required for the solid-state transmitter could 
affect the reliability of the variable capacitors. Problems have been 
encountered on other programs at TRW. When the variable capacitors 
were adjusted too frequently, the movable parts became worn and loose 
and had to be replaced. It was necessary to limit the number of adjust- 
ments which could be made to the variable capacitors. If this transmitter 
is expected to have a long repetitive tuning procedure, steps should be 
taken to prevent this problem with the variable capacitors. 
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Table B-4. Solid State Transmitter 


Part Type 

Quantity 


"9 

Total Stressed Failure Hate (Failures/10 
at 50 Baseplate Temperature 

Hours ) 

Low Capacitor 
Low Transistor 

Low Capacitor 
High Transistor 

High Capacitor 
Low Transistor 

High Capacitor 
High Transistor 

A1 Amplifier stage 

1 





Transistor, power (PT6680) 

1 

489 

754 

489 

754 

Capacitor, variable 

5 

400 

400 

2, 190 

2, 190 

Capacitor, ceramic 

2 

48 

48 

48 

48 

Inductor 

4 

60 

60 

60 

60 

A2 Amplifier stage 

2 





Transistor, power (2N5177) 

1 

638 

983 

638 

983 

Capacitor, variable 

5 

400 

400 

2, 190 

2, 190 

Capacitor, ceramic 

2 

48 

48 

48 

48 

Inductor 

4 

60 

60 

60 

60 

A3 Amplifier stage 

2 





Transistor, power (2N5177) 

2 

1, 504 

2, 318 

1, 504 

2, 318 

Capacitor, variable 

10 

800 

800 

4, 380 

4, 380 

Capacitor, ceramic 

4 

96 

96 

96 

96 

Inductor 

8 

120 

120 

120 

120 

A4 Amplifier stage 

2 





Transistor, power (2N5178) 

3 

2, 598 

4, 004 

2, 598 

4, 004 

Capacitor, variable 

15 

1, 200 

1, 200 

6, 570 

6, 570 

Capacitor, ceramic 

6 

144 

144 

144 

144 

Inductor 

12 

180 

180 

180 

180 

X 3 Multiplier 

2 





Diode, varactor 

2 

1, 270 

1, 270 

1. 270 

1, 270 

Capacitor, variable 

5 

400 

400 

2, 190 

2, 190 

Inductor 

4 

60 

60 

60 

60 

Hybrid {2: 1 or 1:2) 

6 

24 

24 

24 

24 

Hybrid (3*1 or 1:3) 

4 

32 

32 

32 

32 

Miscellaneous resistors, 

20 

40 

40 

40 

40 

metal film 






Total 


20, 129 

25, 524 

46, 979 

52, 374 


NOTES: 


The electrical stresses on the semiconductors were calculated using the power levels given in Figure 2 
A principal electrical stress of 25 percent was assumed for all other parts 



The specifications from Eimac for the 50 -watt TWT include a 
design life of 20, 000 hours. Therefore, for this prediction (28, 600 to 
32, 400 hours MTBF for TWT and converter), it was assumed that the 
TWT would not be operated continuously for more than 20, 000 hours 
during a mission. 
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APPENDIX C 


OSCILLATOR ERROR STATISTICS FOR NAVSTAR 
Introduction 

This appendix summarizes the available information on very stable 
oscillators pertaining to the question of oscillator error in what may be 
called the “autonomous 11 mode of NAVSTAR operation* 

In this mode, the satellite transmissions originate at the satellite 
independently or autonomously by reference to very stable oscillators on 
the satellites as opposed to retransmission of signals from the ground. 
The satellite signals are received at a number of fixed ground stations 
at known locations whose observations are then used to estimate a “total 
ephemeris" consisting of satellite position and velocity elements and 
satellite oscillator bias and one or more oscillator bias rate terms. 

This total ephemeris is then broadcast to all users to allow them to 
correct for differences in satellite oscillators. The bias corrected data 
is then equivalent to fully synchronized satellite transmission. 

To determine optimum bias estimation time constants, ephemeris 
update periods, and the net effect on user position error it_is necessary 
to characterize the time serial statistics of the random oscillator drift. 

Practical considerations of satellite hardware weight, simplicity, 
and reliability have led to an emphasis on the crystal and rubidium 
oscillators so only these two types- are considered here. 

Oscillator Error Characterization 

Consider an oscillator (or clock) providing a phase (or indicated 

time) 

<£{t) 

Relative to a given reference oscillator of frequency co^ (which defines 

co ) we can measure the phase error 
o' 

0*(t) = 0(t) -w o t t 1 ) 


C-l 



f (t) maybe characterized as a random time varying phase error. In 
order to normalize things it is convenient to define the indicated time 
error 


❖ 

t (t) 


/ (t) 


cu 

o 


( 2 ) 


and the associated range error in a system using this as a range delay 
measuring reference 

p*(t) = ct* (t) (3) 

o 

where c is the velocity of light ~ 10 ft/sec. 

Several relatable measures of the essential statistics of these 
three equivalent random functions are in common use. 

1 . Short Time Frequency Stability 

The change in phase error in a time r, divided by t is a 
measure of the frequency error over t 3 thus 

j = £ (t - T) _-_£jt)_ 

t T 


If one measures {t) in N contiguous intervals 


% 

w (t ) 
T O 


W T + T) 




forms the sample mean 


C-Z 



N— 1 


*N <V 




i=0 


= "Nt(‘o + ,N -‘> T ) 


forms the residuals relative to that mean 


Aw (t + It, N) = « r (t + ir) — «' (t ) 
V O ' T ' o ' N'o' 


and the rms residual 


S(t,N) = 


N— 1 


-, 1/2 


N 


■=rj S Aw * 2 ^ + iT ^ 


i=0 


(4) 


Then S (t,N) is the commonly measured short time stability of the 
oscillator . 

2 . Spectrum 

Considered as a random function of time 0 (t) (or equiva- 
lently t or p ) may be characterized in terms of its power spectral 
density function 


<V £ > 

or its autocorrelation function 

V T) 

which by the Wiener-Khintchine theorem is just the Fourier transform 
of G^(t). Usually these are both singular forms since in actuality the 
zero values of both of these functions are infinite. If one is careful 
however, this does not necessarily interfere with their formal usage in 
analysis since the expressions one usually wishes to evaluate m terms 
of G or R are usually finite , i.e., independent in an appropriate 
sense of the zero frequency values of G or R. 
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3. 


Variate Differences 


Given a series of values of <i> at intervals r 

f (t o) 

(t Q + T) 


/(t o + ir) 


form the successive differences 


Aj'*' ^ (t + i-r) = <f>* {t + ir) 
0 o o 


^>,t o + i T , - W ,t 0 ♦ IT, - A «‘) 


etc . 


- + (i - 1)^ ^ 

(t 0 + (i- 1 )^) 


(5) 


and the means and variances 


= <C[ A “ >(t o + iT > " 

Depending somewhat on the type of oscillator the variance of «T itself 
( \Z (i\Z ,n> v9 

(i.e,. , a ^°' } and cr^ 1 are generally indeterminate but cr 

higher order differences generally exist and are finite. 




Relations Between the Various Statistics 

Because of the singularity of some of the functions otherwise in- 
volved, a useful place to start is in terms of the "structure functions," 
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U (t) defined as the mean square difference of phase (or time, or range) 
samples r apart 


U (t) = L* (t + t) - 4* (t)J 

This can be directly related to the autocovariance function, R (t) if the 
latter exists through expanding 

2 2 
U(t) = (t + t)^> + (t)^> 

- 2 < «T(t) «T(t + t)> 

= 2 [r (0) - R(t)] 

R(t) = R (0) ~ 

Further it can be shown (Ref. 2), from the preceeding definitions that 

S 2 ( T ,N) = - 2- |nU(t) - 1u(Nt) 

(N-1)t L N J 

and 

A (2) (t) = 6R(0) - 8R(t) + 2R(2t) 

= 4U (t) - U(2t) 

2 

A (3) (t) = 2 OR (0) - 30R (t ) + i2R (2 t) - 2R (3t) 

= 15U (t) - 6U(2t) + U (3t) etc. 

and finally the phase spectrum is given (if it exists): 
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-jT 


(f ) =4 I R(t) cos 2TrfTdr 


-r 


= 2R (0) 5(f) — 2 J U (t) cos 2-irfTdT 

J 

o 


and ignoring the D.C. term we may write 


r 


(f ) = -2 j U (t) cos Zirirdr (f * 0) 


General Power L>aw Models 

It turns out both experimentally and theoretically that over wide 
ranges of frequency the various functions above may be represented as 
simple powers of their arguments. Thus if we set 

U (r) = Kt V 

we have from the above relations 
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<V f) = 


K_rr 

. vtt r> \ v + * f v+1 - 
r(-v) cos — 2 ~ ^ir) ± 


TTV 


2K r(l + v) sin -j- 

„ iV+ 1 ,v+± _ 
(2-ir) f 


Particular Power Laws 

Particular cases of interest are 

v = -1 corresponding to phase modulation by white noise 

v — i corresponding to frequency modulation by white noise 

v = 2 corresponding to frequency modulation by 1/for 

"flicker" noise 

The first case (v = -1) leads to indeterminate relations between the 
.various measures, and in fact, only G is defined in this case since the 
mean square phase jitter or phase jitter differences of any order are 
infinite. In practice this indeterminacy is resolved by finite circuit 
bandwidths which circumvent the implied ultraviolet catastrophy {Ref. 1). 
However, the details of this resolution need not concern us here since 
over the time constants and frequencies of interest in the NAVSAT prob- 
lem (i.e. , r >> 1 sec) this type of noise is generally not significant. 

The second case (v = 1), corresponding to white noise frequency 
modulation, gives (v = 1) 


s 2 (t,N) 

= {independent of N) 

A (2 > 2 <t> 

= 2Kt 

2 


A< 3 > (t) 

= 6Kt 

<v f > 

2K 



C-7 



The third, case (v = 2) corresponds to flicker frequency modulation 
and is of particular interest since it seems that all oscillators, with the 
possible exception of the maser class, approach such a characteristic 
asymptotically for long times or low frequencies. All the general power 
law relationships above unfortunately go to zero. However, this does 
not correspond to any real physical blow-up, and in fact the ratio of the 
various functions approaches a constant as v approaches 2, which ratio 
may be evaluated by L'Hospitals rule. Normalizing with respect to, say 
arbitrarily, the term 4 — 2 V we have in the limit (v — * 2) 



G ^ f > 


4K 

3~~ 3 

(2ir) f 


sin 


irv 


4 — 2 


n 


Kir 

(2tt) 3 (2 In 2)f 3 


These numerical relationships are summarized below. 


V 

<V f > 

<y f> 

A (2 > 2 (t) 

X(T> 

S 2 (t,N) 

1 

2K 

2K 

2Kt 

6Kt 

K 

(2it ) 2 ? 

T 

2 

K 

K 

Kt 2 

2.43Kt 2 

K /NlnN\ 

109 r 

2 . 78 f 

2.78 yN — ly 
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Comparison of Actual Data 


Figure C-l is a composite plot of several short time frequency 
stability plots for several crystal and one rubidium oscillator from the 
sources indicated. Based on such data and also on McCoubrey's {Ref. 5) 
characterization, we propose to approximate these diverse ensembles by 
the generic representations of Figure C-2. These can be easily trans- 
formed into the spectral domain by the relations previously derived and 

2 9 2 

lead to the effective velocity (proportional through c = (0.983 • 10 ) 
to relative frequency) error spectra as shown on Figure C-3. The quartz 
oscillators at relatively high frequencies show the up-turn spectrum 
effect of white noise phase modulation, introduced through relations de- 
rived in Reference 1 and assuming a respective 160 Hz noise bandwidth; 
but this is really not essential to the present discussion since it only 
affects frequencies above about 1 Hz. 

At low frequencies both quartz and rubidium show the character- 
istic 1/f or flicker frequency or velocity spectra, the rubidium having 
a very slight advantage for frequencies below about 0.001 Hz (with one 
exception to be discussed). 

A major difficulty of the usual short-time stability function at very 
long times concerns the methods of eliminating trends and the dependence 
on N. Usually insufficient attention is paid to these details and the re- 
ported results are somewhat deteriorated for lack of reporting of these 
necessary details . Recognizing these limitations, Barnes (Ref. 3) has 
proposed and utilized the short term variate difference —variance as a 
better behaved statistic. Figure C-4 shows the rms or standard devia- 
tion of variate differences of various orders for a cesium and a quartz 

oscillator. Over the range or values plotted cr varies as r for the 

-3 

quartz oscillator implying a phase jitter spectrum of the form f or 
"flicker" noise in frequency. 

Again the previously derived results permit inferring the frequency 
or velocity error spectra from the preceding results and this is also 
plotted on Figure C-3. It is seen that this oscillator is nearly an order 
of magnitude better than the proposed (Ref. 6) generic characterization 
for quartz oscillators. Several factors may be brought forth to account 
for this difference. 
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Figure C-l. Actual Oscillator Short Time Stability Functions 
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Figure C-3. 


Velocity Error Spectra Due to 
Reference Oscillator Error 
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Figure C-4. Barnes' Variate Difference Data 
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NOT REPRODUCIBLE 









1) This is a Sulzer 5A oscillator which according to Barnes 
at NBS is the best they have ever found. 

2) It has been suggested that as a favored customer NBS may 
well get the cream of the Sulzer production crop. 

3) Constant aging rate effects do not effect the variate differ- 
ence data whereas, depending on the method of removal of 
the trend (unspecified), they may or may not enter the short 
term stability data in Figure C-l on which the generic char- 
acterization is based. 

In verbal conversation with Barnes at NBS it was indicated that the 
best rubidium oscillators show variate differences very similar to those 
of quartz for times of the order of an hour or more, better by a factor of 
perhaps 2, while below t = 1 hour the quartz oscillator is definitely 
superior. 

Direct Interpretation of the Variate Differences 

There is an important direct interpretation of the variate differ- 
ences as given on Figure C-4. Suppose that one performs the following 
expe riment. 

1) Measure n values of oscillator phase, or let us say clock 
time error T, at times 

t = 0, -t, -2t , . . . -{n — 1)t 

2) Fit a polynomial of degree (n— 1) to those n points. 

3) Use the polynomial to extrapolate or predict the time error 
at a time r, one step in the future 

Note that this is quite closely parallel to the problem of calibrating the 
satellite oscillator. The fitting and extrapolation problem can be con- 
veniently expressed in terms of the BaGrange interpolation coefficients, 
thus the extrapolation for the value at 0 in terms of the values at -r , 
-2 t, . . . -nr is given by (Ref. 4, Section 25.2) 


T(0) 


n 


= £ 

K=1 


A^T(-Kt) 
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where 


a: 


n 

K 


= -(-i) 


K 



K = 1, . . . n 


Then the error in extrapolation is 

AT = 'T (0) — T (0) 

But, extending the definition of to include the case K = 0 for which 

A = - 1 we can write 
o 

AT = £ A^T(-Kr) 

K=0 


But this is also exactly the concise formula for the nth difference of the 
T data including the 0th point. (Ref. 4, Section 25. 1) Thus, the rms 
error in extrapolating one t unit forward based on n equally (r) spaced 
points prior is exactly equal to the rms nth difference of the data at in- 
terval t and vice versa. 

Thus, for example, we can read Figure C-4 as saying that for a 
good crystal oscillator, and under equivalent environmental conditions, 
based on a two -day prior (n = 3 point) calibration, we can extrapolate 
time for one day forward with an rms error of 130 nanoseconds or 
130 feet. The second difference (n = 2) data cannot be interpreted in this 
way because this represents not the total rms but the standard deviation 
about a significantly non-zero mean, to determine which requires one 
additional point’ in calibration. 

The reason for these results being so much better than previously 
assumed (conservatively as was then pointed out) is because of the effect 
of aging rate or polynomial component calibration which is not normally 
taken into account but is taken into account exactly in the variate differ- 
ence method.' 
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Simulation 


The important conclusions regarding simulation are summarized 

in Figure C-3 which gives the equivalent range rate spectrum due to 

oscillator frequency error. The frequency range of interest extends 

roughly from the reciprocal sampling time, to the reciprocal data esti- 

2 -fa 

mation time, possibly 10 to 10 Hz. Over this range the quartz error 
spectrum is very nearly 1/f in range rate, while the rubidium spectrum 
is white (f°) at the higher frequencies going into 1/f (note this is power) 
between about 10" to 10" Hz. The white part of this spectrum is easily 
simulated by known methods but no way of simulating 1/f noise in a 
finite state model is presently known. 


References 


1. Mallinckrodt, A.J., "Frequency Synchronization for Doppler 
Tracking Systems," Parti: Stable Oscillators, CRL 498, 

20 April 1965; Part II: Propagation Modes, CRL 509, February 

1966 . 

'2. Allan, D.W., "Statistics of Atomic Frequency Standards," Proc, 
IEEE, Vol. 54, No. 2, p. 221, February 1966. 

3. Barnes, J.A. , "Atomic Timekeeping and the Statistics of Precision 
Signal Aeverators," Proc. IEEE, Vol. 54, No. 2, p. 207, February 

1966. 

4. "Handbook of Mathematical Functions," AMS 55, National Bureau 
of Standards, U.S. Government Printing Office. 

5. McCoubrey, A. O. , "A Survey of Atomic Frequency Standards," 
Proc. IEEE , Yol. 54, No. 2, p. 116, February 1966. 

6. "Study of Navigation and Traffic Control Employing Satellites," 
Proposal No. 8710 . 000 to NASA-ERC , TRW Systems, January 24, 

1967. 


C- 16 



